Hyperbolic Polynomials Approach to Van der 
Waerden/Schrijver- Valiant like Conjectures : 
Sharper Bounds , Simpler Proofs and Algorithmic 

Applications 

Leonid Gurvits * 
February 2, 2008 

Abstract 

Let p{xi, ...,x„) = p(X),X £ i?" be a homogeneous polynomial of degree n in n real 
variables , e = (1,1,..,1) G i?" be a vector of all ones . Such polynomial p is called e- 
hyperbolic if for all real vectors X e i?" the univariate polynomial equation p{te — X) — 
has all real roots Ai(X) > ... > Xn{X) . The number of nonzero roots |{i : Xi{X) ^ 0}| is 
called Rankp{X) . A e-hyperbolic polynomial p is called POS'-hyperbolic if roots of vectors 
X S i?" with nonnegative coordinates are also nonnegative (the orthant i?" belongs to the 
hyperbolic cone) and p{e) > . Below {ei, e„} stands for the canonical orthogonal basis 
in R". 

The main results states that if p{xi,X2, x„) is a POS'-hyperbolic (homogeneous) polyno- 
mial of degree n , Rankp{ei) — Ri and p{xi, X2, Xn) > ni<i<n a;^; Xi > 0, 1 < i < n, 
then the following inequality holds 

^^^P(0,...,0)> n (^)--(G.=n.in(P.n+l-.)). 

l<i<n 

This theorem is a vast (and unifying) generalization of the van der Waerden conjecture 
on the permanents of doubly stochastic matrices as well as the Schrijver- Valiant conjecture 
on the number of perfect matchings in /c-regular bipartite graphs . These two famous results 
correspond to the POS'-hyperbolic polynomials being products of linear forms. 

Our proof is relatively simple and " noncomputational" ; it actually slightly improves 
Schrijver's lower bound , and uses very basic ( more or less centered around RoUe's theorem 
) properties of hyperbolic polynomials . 

We present some important algorithmic applications of the result, including a polynomial 
time deterministic algorithm approximating the permanent of n x n nonnegative entry-wise 
matrices within a multiplicative factor ^ for any fixed positive m . This paper introduces 
a new powerful "polynomial" technique , which allows as to simplify /unify famous and hard 
known results as well to prove new important theorems . 

The paper is (almost) entirely self-contained , most of the proofs can be found in the 
Appendices. 
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1 Introduction 



An n X n matrix A is called doubly stochastic if it is nonnegative entry-wise and every column 
and row sum to one. The set of n x n doubly stochastic matrices is denoted by Q^i- Let 
A(A;, n) denote the set of n x n matrices with nonnegative integer entries and row and column 
sums all equal to k . We define the following subset of rational doubly stochastic matrices : 
^k,n = {k-^A : A G A(fe,n)} . In a 1989 paper .5^ R.B. Bapat defined the set Dn of doubly 
stochastic n-tuples of n x n matrices. 

An n-tuple A = (^i, • • • , An) belongs to Dn iff Ai y 0, i.e. Ai is a positive semi-definite matrix, 
1 < i < n ; trAi = 1 for 1 < i < n ; J2i=i Ai = I, where /, as usual, stands for the identity 
matrix. Recall that the permanent of a square matrix A is defined by 

n 

per{A) = l[A{i,a{i)). 

Let us consider an n-tuple A = (^i, A2, ...An), where Ai = {Ai{k, I) : 1 < k,l < n) is a complex 
n X n matrix (1 < i < n). Then detQ2i<i<n ^iAi) is a homogeneous polynomial of degree n in 
ti,t2, - ■ ■ , tn. The number 

M(A) := D{Ai,A2, ■ ■ ■ , A^) = ^ det(ti^i + • • • + (1) 

Oti--- dtn 

is called the mixed discriminant Ai, A2, ■ ■ ■ , An . 

The permanent is a particular (diagonal) case of the mixed discriminant . I.e. define a multilin- 
ear polynomial Mn/ a (ti, tn) = ni<i<nEi<j<n^(*'i)*i- Thenper(^) = Q^^^MulA{h, ...,tn). 
Let us recall two famous results and one recent result by the author. 

1. Van der Waerden Conjecture 

The famous Van der Waerden Conjecture jlj states that minAennD^A) = ^ (VDW- 
bound) and the minimum is attained uniquely at the matrix J„ in which every entry 
equals ^. Van der Waerden Conjecture was posed in 1926 and proved only in 1981 : D.I. 
Falikman proved in [13] the lower bound ^ ; the full conjecture , i.e. the uniqueness part 
, was proved by CP. Egorychev in ^21 ■ 

2. Schrijver- Valiant Conjecture 

Define 

X{k,n) = min{per(^) : A G 9.k,n} = fc"" min{per(A) : A € Ak,n};0ik) = lmi^{X{k,n))n . 

It was proved in j^l] that , using our notations , 9{k) < g{k) = (^^)'^~^ and conjectured 
that 9(k) = g{k) . Though the case of fc = 3 was proved by M. Voorhoeve in 1979 [36.. , 
this conjecture was settled only in 1998 |35j (17 years after the published proof of the Van 
der Waerden Conjecture). The main result of is the following remarkable inequality : 
min{per(^) : A G J7fc,„} > (^)('=-i)" (Schrijver-bound) . 

The proof in is probably one of the most complicated and least understood in the 
theory of graphs . 
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3. Bapat's Conjecture (Van der Waerden Conjecture for mixed discriminants) 

One of the problems posed in [H] is to determine the minimum of mixed discriminants of 
doubly stochastic tuples : minA&D„D{A) =? 

Quite naturally, R.V.Bapat conjectured that minA£D„D{A) = ^ (Bapat-bound) and 
that it is attained uniquely at J„ =: (^1, ^/). 

In |51 this conjecture was formulated for real matrices. The author had proved it |HJ 
for the complex case, i.e. when matrices Ai above are complex positive semidefinite and, 
thus, hermitian. 

The (VDW-bound) is the simplest and most powerful bound on permanents and therefore 
among the simplest and most powerful general purpose bounds in combinatorics. Besides its 
many applications to the graph theory and combinatorics , (VDW-bound) has been recently 
used for deterministic approximations of permanents .(Much more recent proof of (Bapat- 
bound) was actually motivated by the scaling algorithm ^H] ,^^1 to approximate mixed dis- 
criminants and mixed volumes.) It is easy to check that (Schrijver-bound) is implied by 
(VDW-bound) for A; > n : ^ = ni<fc<n(^)''"^ > (^)('=-i)" . Therefore , it was not clear 
whether the scaling algorithm in gives better approximating exponent for sparse matrices 
: the "scaled" doubly stochastic matrix may have irrational entries even if the input matrix is 
boolean and (Schrijver-bound) is superior to the (VDW-bound) only on "very" rational 
sparse doubly stochastic matrices . 

Since our generalized (Schrijver-bounds) (18), (19) depend only on the "sparsity" hence the 
scaling algorithm for permanents in indeed gives better approximating exponent for sparse 
matrices (scaling algorithm for mixed discriminants in JH] ,^^1 gives better approximating 
exponent for tuples of "small" rank PSD matrices). 

1.1 Van der Waerden / Schrijver- Valiant like conjectures and homogeneous 
polynomials 

Let Hom{m,n) be the linear space of homogeneous polynomials p{x),x G EJ^ of degree n 
in m real varibles ; correspondingly Homj^{ra^n){Horaj^j^{rn,n)) be a subset of homogeneous 
polynomials p{x),x G i?™" of degree n in m real varibles and nonnegative(positive) coefficients 

Definition 1.1: 

1. Letp G Hom+(n,n),p{xi,...,Xn) = E(ri,...,r„)6/n,„ «(r-i,...,r„) ni<i<n 4' be a homogcneous 
polynomial of degree n in n real variables. Here Im,n stands for the set of vectors r = 
(ri, ...,rm) with nonnegative integer components and X]i<i<m^« = 

The support of the polynomial p{xi, ■■■,Xn) as above is defined as supp{p) = {(ri, r„) G 
n '■ 0'{ri,...,rn) 7^ 0} . The convex hull CO{supp[p)) of supp{p) is called the Newton 
polytope of p . 

For a subset A C {l,...,n} we define Sp{A) = max(r^^ j,^)g^„pp(p)X]jgA'^«- Given a 
vector (oi, a^) with positive real coordinates , consider univariate polynomials Dji(t) = 
p{t{J2ieAei) + J2i<j<n'^jej),VA{t) = p{t{J2ieAei) + J2jeA' aj^j) • Then S'p(^) is equal to 
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the degree of the polynomials Da, VA{t) : 

Sp{A) = deg{DA) = deg{VA) 



(2) 



2. The following linear differential operator maps Hom{n,n) onto Hom(n — l,n — 1) : 

d 

We define Pxjj 2 < z < n in the same way for all polynomials p G Hom{n, n). Notice that 

p{xi, ...,Xn) = XiPa:^{x2, ■.■,Xn) + q{xi, ...,Xn);qxi = ^- (3) 

The following inequality follows straight from the definition : 

Sp^^{A) < mm{n - 1, Sp{A)) : A C {2,...,n},p G Hom+{n,n). (4) 

3. Consider p G Homj^{n,n) We define the Capacity as 

Cap{p)= inf p{xi,...,Xn). 

2^'>0'ni<i<„^i=i 

It follows that if p G Hom-\-{n, n) then 
Notice that 

log{Cap{p)) = infY y^=o log(p(e^% e^")), 
and if p G Hom+{n, n) then the functional log(p(e^i , e^")) is convex . 

4. Consider a stratified set of homogeneous polynomials : F = Ui<n<oo , where Fn G 
Hom-\.{n,n) . We call such set VDW- FAMILY if it satisfies the following properties : 

(a) If a polynomial p G F„, n > 1 then for all 1 < z < n the polynomials p^;. G Fn-i- 
(b) 

Cap{px,) > g{Sp{{i]))Cap{p) ■.p€Fj,l<i< j; g{k) = k>l. (6) 

I 

Example 1.2: Let A = {A{i,j) : 1 < i < n} he n x n matrix with nonncgativc entries . As- 
sume that ^i<j<n^ihj) > for all 1 < i < n. Define the following homogeneous polynomial 
MulA{ti,...,tn) = Ui<i<nT,i<j<nMiJ)tj ■ Clearly , MuIa G Hom+{n,n) and MuIa / . 
It is easy to check that SmuIaHj}) = I{^ • ^ihj) 0}| {SMulAiij}) equal to the number of 
non-zero entries in the jth column of A) . 

Notice that if A G A{k,n) (or A G i}.{k,n)) then SmuIaHj}) <k,l<j<n . 

More generally , consider a ra-tuple A = (Ai,A2, ...An) , where the complex hermitian n x n 
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matrices are positive semidefinite and J2i<i<n ^ (their sum is positive definite). Then the 
homogeneous polynomial DETA{ti, ■■■,tn) = det(X]i<i<n ^i^i) ^ Hom+{n,n) and DETa / . 
Similarly to polynomials MuIa , we get that 5'd£;ta({j}) = Rank{Aj), I < j < n. 
The Van Der Waerden conjecture on permanents as well as Bapat's conjecture on mixed discrim- 
inants can be equivalently stated in the following way (notice the absence of doubly stochasticity 
): 

The van der Waerden conjecture on the permanents corresponds to polynomials MuIa G 
Hom+{n,n) : A > , the Bapat's conjecture on mixed discriminants corresponds to DETa G 
Hom^{n,n) : A ^ . The connection between inequality (7) and the standard forms of the 
van der Waerden and Bapat's conjectures is established with the help of the scaling (^7] , ^H] 
, ^ni)- Notice that the functional log(p(e^i, ...,e^^)) is convex if p G Hom^{n,n). Thus the 
inequality (7) allows a convex relaxation of the permanent of nonnegative matrices and the 
mixed discriminant of semidefinite tuples . This observation was implicit in ^Jj and crucial in 

I 

1.2 The Main (polynomial) Idea 

The following (meta)theorem describes the main idea of this paper . 

Theorem 1.3: Let F = Ui<n<oo -^n ^ VDW-FAMILY and the homogeneous polynomial 
p & Fn. Then the following inequality holds : 

n ff(min(5p({i})),n + 1 - i))Cap{p) < -— p(0, ...,0) < Cap{p). (8) 

l<i<n " 

Proof: Our proof is by natural induction . Notice that the function g{k) = (^j^)^~^ is strictly 
decreasing on the semiinterval [1,cxd) (we define ^(0) = 1). The theorem is obviously true 
for n = 1 . Suppose it is true for all n < A; < cxo and the polynomial p £ Fk+i ■ Since 
F = Ui<n<oo Fn is a VDW-FAMILY hence p^,^, G F^; Cap{p,,^,) > g{Sp{{k + l}))Cap{p) . 
It follows from (obvious) inequality (4) that Sp^^^^{{i}) < Sp{{i}), 1 < i < k . Therefore , by 
induction , we get the needed inequality : 

CT^^(O'---'O) = asfT7(5^I^P-.+i(0'---'0)) > 

> ni<i<feS'(min(S'p^^^^({i})),n + 1 - i))Cop(p^,^J > 

> ni<^<fe g{mm{Sp^^^^ {{i})),n + 1 - i))Cap{p)g{Sp{{k + 1})) > 

> ni<i<fc+i 5'(min(5p({i})), n + 1 - i))Capip). | 

Corollary 1.4: 
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1. If the homogeneous polynomial p ^ then 

2. If the homogeneous polynomial p G Fn and Sp{{i})) < k < n,l < i < n then 

I, _ 1 LI an 

( )(fc-i){n-fc)_Cap(|5) < —p{0,...,0) < Cap{p). (10) 

K k'^ OXl...OXn 

Proof: Both inequalities follow the main inequality (8) and from the next easily proved identity 

^= n 9{k). (11) 

l<fc<n 

I 

What is left now is to present a VDW-FAMILY which contains all polynomials DETx 
, where the n-tuple A = {Ai, An) consists of positive semidefinite hermitian matrices (and 
thus contains all polynomials MuIa , where ^ is n x n matrix with nonnegative entries). If such 
VDW-FAMILY set exists then the Van der Waerden , Bapat , Schrijver- Valiant conjectures 
would follow (without any extra work , see Example 1.2) from Theorem 1.3 and Corollary 1.4 . 
One of such VDW-FAMILY , consisting of POS'-hyperbolic polynomials , is defined in the 
next section. 

2 Hyperbolic polynomials 

The following concept of hyperbolic polynomials was originated in the theory of partial differ- 
ential equations llSf . jIBI ■ recently became "popular" in the optimization literature 111^ 
. \1CII .]S^ ^ . The paper gives nice and concise introduction to the area (with much simplified 
proofs of the key theorems ) . 

Definition 2.1: 

1. A homogeneous polynomial p : C" ^ C of degree n{p € Hom{m, n)) is called hyperbolic 
in the direction e E EJ^ (or e- hyperbolic) if p{e) ^ and for each vector X € R"^ the 
univariate (in A) polynomial p{X — Ae) has exactly n real roots counting their multiplic- 
ities. 

2. Denote an ordered vector of roots oip{x — Xe) as X{X) = (An(V) > An-i(V) > ...Ai(V)). 
Call X £ e-positive (e- nonnegative) if Ai(V) > (A„(V) > 0). We denote the closed 
set of e- nonnegative vectors as Ne{p), and the open set of e-positive vectors as Ce{p)- 
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Definition 2.2: Let p : C"^ C he a homogeneous polynomial of degree n in m variables. 
Following [2ni , we define the p- mixed form of an n- vector tuple X = (Xi, .., X„) : Xi G C" as 

Mp(X) =: Mp{Xi, ..,Xn) = Q^^^^Pi E (12) 
The following polarization identity is well known 

Mp{Xi,..,Xn) = 2-'' J2 P(T. ^i^i) n (13) 

6ie{-l,l},l<j<n l<i<n l<i<n 

Associate with any vector r = (ri, ...,r„) G In.n an n-tuple of m-dimensional vectors con- 
sisting of ri copies of < z < n). It follows from the Taylor's formula that 

p{ y: = E n «?M,(x,) ^ (14) 

l<i<n rg/„,„ l<j<n lll<j<n'«- 

I 

We collected in the following proposition the properties of hyperbolic polynomials used in this 
paper . 

Proposition 2.3: FACT 1 . 

p{X)=p{e) n (15) 

l<i<n 

FACT 2 . If p is e-hyperbolic polynomial and p{e) is a real nonzero number then the coeffi- 
cients of p are real (J^ , follows from (15) via the standard interpolation). If p is e-hyperbolic 
polynomial and p{e) > then p{X) > for all e-positive vectors X G Ce{p) C i?"* . 

FACT 3 . Let p £ Hom{m,n) he e-hyperbolic polynomial and d G Ce{p) C . Then p 
is also d- hyperbolic and C^ip) = Ce{p), N^ip) = Ne{p) . fllS}/ , \2^^ , very simple proof in TSSj .) 

FACT 4 . Let p G IIom(m,n) be e-hyperbolic polynomial . Then the polynomial Pe{X) =: 
^p{X + te)|(f=o);Pe G IIom{m,n — 1) is also e- hyperbolic and Ce{p) C Ce{pe) (W^ > W^i , 
Rolle's theorem ). 

FACT 5 . Let p G Hom{m,n). Then the p-mixed form AIp{Xi, .., Xn) is linear in 
each vector argument Xi G C"^. Let p G IIom(m,n) be e- hyperbolic and p{e) > . Then 
Mp{Xi, ..,Xn) > if the vectors Xi G R"^, 1 < i < n are e-positive , proved by induction 
using FACT 4) . 

We use in this paper the following sub-class of hyperbolic polynomials . 

Definition 2.4: A polynomial p G Hom{m,n) is called P05-hyperbolic if p{e) > 0, e = 
(1,1,...,1) G R"^ ; p is e-hyperbolic and the closed convex cone Ne{p) contains the nonnegative 
orthant i?™ . (In other words , all the roots of the univariate polynomial equation p{X — te) = 
are nonnegative if the coordinates of the vector X are nonnegative real numbers.) 
It follows from the identity (14) and FACT 5 that P05-hyperbolic polynomials have nonneg- 
ative coefficients . I 
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Probably the best known example of a hyperbolic polynomial comes from the hyperbolic 
geometry : p{xQ, ...,Xk) = Xq — Yl,i<i<k^i- This polynomial is hyperbolic in the direction 
(1, 0, 0, 0). Another "popular" hyperbolic polynomial is det(X) restricted on a linear real 
space of hermitian n x n matrices . In this case mixed forms are just mixed discriminants 
, hyperbolic direction is the identity matrix I , the corresponding closed hyperbolic cone of 
/-nonnegative vectors coincides with a closed convex cone of positive semidefinite matrices . 
Less known , but very interesting , hyperbolic polynomial is the Moore determinant Mdet(y) 
restricted on a linear real space of hermitian quaternionic n x n matrices . (The Moore deter- 
minant is a particular case of the generic norms on Jordan Algebras.) The Moore determinant 
is , essentially , the Pfaffian (see the corresponding definitions and the theory in a very readable 
paper ^40 ) . 

This paper benefits from the fact that as multiUnear polynomials MuIa G Homj^{n, n) : 
j4 > 0, > , as well determinantal polynomials DETa G Homj^{n,n) : A ^ 
0, Ei<i<n ^ are POS-hyperbolic. 

2.1 FOS'-Hyperbolic polynomials form VDW- FAMILY 

Let q E Hom+{n,n) be a P05-hyperbolic polynomial. For a vector X G C" we define the 
integer number Rankq{X) as the number of nonzero roots of the equation q{X — te) = 0, e = 
(1,1,..., 1) = Z]i<i<n^«- It follows from the identity (2) that Rankq{J2i£A^i) — Sq{A),A C 
{1,2,. 

Theorem 2.5: 

1. Let q G Hom^{n,n) be PO S -hyperbolic polynomial . If 1 < Rankq{ei) = k <n then 

Cap{q^,)>g{k)Cap{q) (16) 

2. Let g(xi, X2, x„) be a POS -hyperbolic (homogeneous) polynomial of degree n . Then 
either the polynomial = or q^^ is POS -hyperbolic . If Cap{q) > then qx^ is 
(nonzero) POS -hyperbolic . 

Corollary 2.6: Let PHP{n) C Hom^{n,n) be a set of POS -hyperbolic polynomials of degree n 
in n variables ; define PHP+{n) = {p G PHP{n) : Cap{p) > 0}. Then as U„>i(Pi?P(n)U{0}) 
as well Un>iPHP+{n) is VDW- FAMILY . 

The second part of Theorem 2.5 is , up to minor modifications , well known (FACT 4; see , 
for instance, [23]) IHH] ) and follows from Rolle's theorem . The main new results "responsible" 
for the first part of Theorem 2.5 are the next Lemma 2.7 and its Corollary 2.8. 

Lemma 2.7: 
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1. Let ci, c„ be real numbers ; < Ci < 1,1 < i < n and J2i<i<n (H = n — 1. 
Define the following symmetric functions : 

l<i<n l<i<nj^i 

Then the following entropic inequality holds : 

2. (Mini van der Waerden conjecture) 

Consider a doubly-stochastic n x n matrix A = [a|6|...|6] . I.e. A has n — 1 columns 
equal to the column vector b , and one column equal to the column vector a . Let a = 
(oi,...,a„)^ : Oi > 0,J2i<i<nO'i = 1 ; b = (6i,...,5„)^ : bi = < i < n. Then the 

permanent Per (A) > ^ . 

Corollary 2.8: Consider an univariate polynomial 

R{t) = X]o<i<rt^«** = Ili<i<ni^i^ + ^i) ) where ai,bi > . If for some positive real number C 
the inequality R{t) > Ct holds for all t > then 

di = ^R{0) > C((^r-i) (17) 
The inequality (17) is attained only on the polynomials R{t) = A(t + a)"' : A,a > 0. 



2.2 Newton Inequalities , Alternative Proof of Corollary 2.8 , Volume Poly- 
nomials 

Let R{t) = X]o<i<n dit^ be an univariate polynomial with real coefficients . If such polynomial 
R has all real roots then its coefficients satisfy the following Newton's inequalities : 

NIs : dj > di-idi+i „ : 1 < i < n - 1. 

The following weak Newton's inequalities WNIs follow from NIs if the coefficients are non- 
negative: 

WNIs ■ didt^ < — I ^ I : 2 < z < n. 
n \i I 

Lemma 2.9: Let R{t) = J2o<i<ndit'^ ^'^ univariate polynomial with real nonnegative co- 
efficients satisfying weak Newton's inequalities WNIs -If for some positive real number C the 
inequality R{t) > Ct holds for allt>0 then 

di > 
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Proof: If do = then di > C > C((^)'^-^) . Thus we can assume that do = 1 . It fohows 
from weak Newton's inequahties WNIs that 



n \il 



2<i<n. 



Therefore for nonnegative values of t > we get the inequahty 



n \1 I n \2 I n 



n n 



Which gives the inequahty (1 + > Ct . The inequality di > C{{^^)'^ ^) follows now 

easily . I 



Remark 2.10: The Newton Inequalities are not sufficient for the real rootedness . The classi- 
cal example is provided by some univariate volume polynomials R{t) = Vol{tCi + C2) , where 
Ci,C2 are convex compact sets (see , for instance , j^H]) • In this case the Newton Inequalities 
follow from the celebrated Alexandrov-Fenchel Inequalities. 

Using Lemma 2.9 , the Alexandrov-Fenchel Inequalities and a bit of extra work allows to 
extend the results of this paper , i.e. Theorem 1.3 ,to the multivariate volume polynomials 
Vol(tiCi + t2C2 + ... + tnCn) , where Ci,C2,...,C„ are convex compact subsets of R"" . In 
other words , there exists a VDW-FAMILY which contains all such volume polynomials 

VolihCi+t2C2 + ...+tnCn) ■ 

This extension leads to a randomized poly-time algorithm to approximate the mixed volume 
M(Ci, Cn) = dti'^"i)t„ ^ ^K'tiCi + ^2^2 + ... + tnCn) within exponential factor e" . The algo- 
rithm is pretty much the same as in Theorem 4.7 in this paper , the randomization is needed 
to evaluate the oracle , i.e. to evaluate Vol{tiCi + t2C'2 + ... + tnCn) ■ 

The best current appoximation factor is n'^^") ([311 , (randomized) ; ^S] ) ^S] (determin- 
istic) . I 



Definition 2.11: Let p G Hom+{n,n),p{xi, ...,Xn) = E(ri,...,rn)G7„.n ni<i<n be 

a homogeneous polynomial with nonnegative coefficients of degree n in n real variables. Call 
such polynomial AF-Polynomial if the following Alexandrov-Fenchel Inequalities hold : 

Mp(Xi, X2, X3, Xn)'^ > Mp{Xi,Xi,X3, ...,Xn)Mp{X2,X2,Xs, Xn) Xi, X„ G i?" . 

(The p-mixed form Mp{Xi, X2, X3, ...,Xn) is defined in Definition 2.2 (formula (12) .) 
We denote as AF(n) a set of all AF-Polynomial of degree n in n real variables and define 
AF^{n) : {p G AF{n) : Cap{p) > 0} . We denote as Vol{n) a set of polynomials Vol{tiCi + 
t2C2 + ... + tnCn) , where Ci, C2, C„ are convex compact subsets of . 

Notice that PHP[n) C AF{n) ( POS -Hyperbolic polynomials are AF-PolynomialsJ 12!^ 
; Vol{n) C AF{n) , this inclusion is just a restatement of the celebrated Alexandrov-Fenchel 
Inequalities for mixed volumes 

I 
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Theorem 2.12: as Un>iAF{n) as well U„>iAF+(n) is VDW-FAMILY . 

Proof: The definition of the VDW-FAMILY consists of two properties (see Part 4 of Defini- 
tion 1.1) . The property (a) follows from the definition of AF-Polynomials : 

d 

Pxi{x2,-,Xn) = ^— p(0,X2,...,Xn) = ((n - 1)1)"^ Mp{ei,Y,Y, ...,¥); Y = {0,X2, ...,Xn). 
ox I 

The property (b) follows from Lemma 2.9 . Indeed , if X,Y £ i?" and p is a AF-Polynomial 
then the coefficients of the univariate polynomial p{tX + Y) are nonnegative and satisfy the 
Newton inequalities NIs . I 



3 Harvest 

Corollary 2.6 allows to "plug-in" POS'-hyperbolic polynomials to Theorem 1.3 . The most 
spectacular application is the following generalization of (Schrijver-bound) . 

Theorem 3.1: Let A = {A{i,j) : I < i,j < n} be a matrix with nonnegative entries . Define 
Cj = Card{{i : A{i,j) ^ 0}), 1 < j < n. I.e. Cj is the number of non-zero entries in the jth 
column of A . Then 

per{A) > Yl 5(min(Cj, n + 1- j))Cap{MulA)- (18) 

l<j<n 

If Cj < k,l < j < n then 

per{A) > (^^)ik-i)i--k)^Cap{MulA). (19) 
(Recall that if A is doubly stochastic then Cap{MulA) = ^ ■) 

Remark 3.2: The lower bound (18) can be viewed as a NONREGULAR generalization 
of (Schrijver-bound) ; it "interpolates" between (VDW-bound) {Cj = n) and the sparse 
case {Cj « n) . The lower bound (19) is actually sharper than (Schrijver-bound) : |^ = 

ni<,<.5(i)><?(fe)^ = (¥)^'=-')('=). ■ 



4 Algorithmic Applications 



Suppose that a POS'-hyperbolic polynomial 

p{xi,...,xn)= X! «(n,.-,rn) n 

has nonnegative integer coefficients and is given as an oracle . I.e. we don't have a list coeffi- 
cients , but can evaluate p{xi, ...,Xn) on rational inputs . 
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A deterministic polynomial-time oracle algorithm is any algorithm which evaluates the given 
polynomial p{.) at a number of rational vectors g*-*^ = {qi \ ■■■,Qn'^) which is polynomial in n 
and log(p(l, 1, .., 1)); these rational vectors g^*-* are required to have bit-wise complexity which 
is polynomial in n and log(p(l, 1, .., 1)) ; and the number of additional auxilary arithmetic com- 
putations is also polynomial in n and log(p(l, 1, .., 1)) . 

If the number of oracle calls ( evaluations of the given polynomial p{.)) , the number of addi- 
tional auxilary arithmetic computations and bit-wise complexity of the rational input vectors 
are all polynomial in n (no dependence on log(p(l, 1, .., 1)) ) then such algorithm is called 
deterministic strongly polynomial-time oracle algorithm. 

The following result was proved in j39j . 
Theorem 4.1: 

1. Letp G Hom^{n, n) he POS -hyperbolic polynomial . Then the function Rankp{J2i£A ~ 
Sp{A) is submodular , i.e. Sp{ALlB) < Sp{A) + Sp{B) - Sp{An B) : A,B C {1,2, ...,n} . 

2. Consider a nonnegative integer vector r = (ri, ■■■,rn),J2i<i<n^i = ''^ • Then r € supp{p) 
iffr{S) = Y.i^sn < Sp{S) : S C {1,2, . 

Corollary 4.2: Letp G Hom+{n,n) be POS -hyperbolic polynomial. Associate with this poly- 
nomial p the following bounded convex polytope : 

SU Bp = {{xi,...,Xn) -.^Xi < Sp{S) : S d {1,2,.. .,n]; ^ = n; > 0, 1 < f < n} 

i<^S l<i<n 

Then SUBp is equal to the Newton polytope of p , i.e. SUBp = CO{supp{p). 

Proof: The inclusion CO{supp{p)) C SUBp follows from the definition of the function Sp . 
Since the function Sp : 2'f^'^''' '"'^ {0, 1,2, ...,n} is submodular and integer valued hence the 
extreme points of the polytope SUBp are integer nonnegative vectors ^1] . Using the second 
part of Theorem 4.1 , we conclude that all the extreme points of the polytope SUBp belong to 
supp{p) . It follows from the Krein-Milman Theorem that SUBp C CO{supp{p)). | 

Corollary 4.3: Given POS -hyperbolic polynomial p € Hom+{n,n) as an oracle , there exists 
strongly polynomial-time oracle algorithm for the membership problem as for supp{p) as well 
for the Newton polytope CO{supp{p)). 

Proof: Let X = (xi, be a vector with real nonnegative coordinates , sumi<i<nXi = n 

. Consider a function Gx{S) = Sp{S) — J2ies^i^^ ^ {l,2,...,n}. Then Gx is submodular 
and X G CO{supp{p)) iff min5(-{i,2,...,n} G'x(5') > (X G supp{p) iff min5(-{i,2,...,n} GxlS) > 
and X is integer). In view on the recent results on the minimization of submodular functions 
we only need to prove that there exists a strogly polynomial-time oracle algorithm to compute 
Gx{S) . Computing X^ieS'^* ^^^^ • ^"^^ Sp{S) = deg{DA) , where the univariate polynomial 
DAit) = pitiJ2ieA^i) + J2i<j<n^j) ■ Clearly we can compute the degree deg{DA) via the 
standard interpolation , which amounts to at most n + 1 evaluations of p and 0(nlog(n)^ 
arithmetic operations . I 
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Remark 4.4: Consider p G i?om+(4,4),p = xiX2X3X4,+xlxl. Then S'p({l, 2, 3}) = 3, 5p({l,2}) = 
2, Sp{{2, 3}) = 2, Sp{{2}) = 2 and therefore the function Sp is not submodular . It is easy to see 
that Sq is submodular for all q G Hom+{3, 3). On the other hand , consider p G Hom+{3, 3),p = 
X1X2 + + X3X1. One can check that SUBp ^ CO{supp{p). The proof in of submod- 
ularity of Sp for P05-hyperbolic polynomials p is based on the proved LAX conjecture , 
0, ISZI)) very nonelementary result . 

It was proved in j39j that , unless P = NP , there is no deterministic polynomial-time oracle 
algorithm to check if (1, 1, 1) G supp{p) for integer polynomials p € Homj^{n,n) . | 



Definition 4.5: A homogeneous polynomial q G Homj^{n,n) is called doubly-stochastic if 
1, 1) = 1 for all 1 < i < n. The doubly-stochastic defect of the polynomial q is 
defined as 

^5(g) = Ei<.<n(a|-9(l'l' •••>!) -I)'- 

A polynomial q G Hom^(n, n) is called scalable if there exists a positive vector (3 = f3n) 
such that the scaled polynomial (^/^(xi, = q{J2i<i<n Pi^i) is doubly-stochastic . (It is 
easy to see that q G Hom+{n,n) is scalable iff the infimum inf^.^Q-rr ^ =ili^^^ ■■•^^n) = 

Cap{q) is attained .) A polynomial q G Hom^{n,n) is called indecomposable if infimum 
inf^^^Qj-j Xi=i 9(^1' ~ Cap{q) is attained and unique . (Theorem D.l in Appendix 

D "justifies" , in the P05-hyperbolic case , our notion of indecomposability.) I 



Theorem 4.6: A POS -hyperbolic polynomial q G Hom^{n,n) is indecomposable if and only 

if the following two equivalent conditions hold : 

Condition 1. ^ ^ . s— 9(0, 0) > : 1< i / j < n . 

Condition 2. Rankp{J2ieA ^i) = Sp{A) > |^| : yl C {1, 2, n}, 1 < |^| < n . 



The following theorem combines the algorithm and its analysis from (see section 4 jl9jMn 
and Theorem . Similarly to jl9j , we use the ellipsoid method to approximate miny^ ^,^q log(p(e 

. The starting ball is centered at and has the radius n2 log(2p(l, 1, 1)) . To run the ellipsoid 
method we need to compute the gradient of log(p(e°^, e"") ; since p is homogeneous polyno- 
mial of degree n hence we can compute the gradient by n standard univariate interpolations. 
These univariate interpolations amount to O(n^) oracle calls and 0{v?) arithmetic operations 
. The ellipsoid updating also requires 0{v?) arithmetic operations . 



Theorem 4.7: There exists a deterministic polynomial-time oracle algorithm which computes 
for given as an oracle indecomposable POS -hyperbolic polynomial p{xi, .■■,Xn) a number F{p) 
satisfying the inequality 

p(0,...,0) <F(p) <2( n 9imm{Spi{i})),n + l-i)))-^- p{0,...,0)< 



dxi...dXn l<i<n dxi...dXn 

~ n\ dxi...dxr. ' ' ^ 



12 



Theorem 4.7 can be (slightly) improved . I.e. it can be applied to the polynomial 



Notice that the polynomial pk is a homogeneous polynomial of degree n — k in n — k variables 
. It follows from Theorem 2.5 that ii p = pq is P05'-hyperbolic and Cap{p) > then for all 
< k < n the polynomials pk are also POS'-hyperbolic and Cap{pk) > . Also , if p = po is 
indecomposable then pk is indecomposable as well (Theorem 4.6). 

The trick is that if k = mlog2{n) then (using the polarizational formula (13) ) the polynomial 
Pk can be evaluated using 0(n™+^) oracle calls of the (original) polynomial p = po . This 
observations allows to decrease the worst case multiplicative factor in Theorem 4.7 from to ^ 
for any fixed m . If the polynomial p = po can be explicitly evaluated in deterministic polynomial 
time , this observation results in deterministic polynomial time algorithms to approximate 
dxf dx 0) within multiplicative factor ^ for any fixed m . Which is an improvement of 

results in jl7j (permanents , p is a multilinear polynomial) and in jl8j , [15] (mixed discriminants, 
p is a determinantal polynomial) . 



Remark 4.8: Let A be an n x n matrix with nonnegative entries , S C {1, 2, n}, l^l = 
m log2(n) . Assume , modulo polynomial time preprocessing , that A is fully indecomposable J7j 
. Using the Laplace expansion for permanents , we get that per{A) = J2\T\=n-\s\ per{As^T)per{As' 
. This suggest the following deterministic algorithm : compute exactly the permanents per{As^T) 
of "small" matrices As^t and run the algorithm from to approximate with the multiplicative 
factor ^ the permanents of "large" matrices As'^t' ■ This algorithm achieves the multiplica- 
tive factor ^ , but it runs in quasi-polynomial time . Our approach is to apply Theorem 
4.7 to indecomposable P05-hyperbolic polynomial J2\T\=n-\s\P^^i^s,T)MulAg, ^, , which can 
be evaluated in deterministic Po/y(n)-time. Our new "hyperbolic" (VDW-bound) (9) al- 
lows multiplicative factor ^ . We can use the same trick for sparse matrices using our new 
"hyperbolic" (Schrijver-bounds) (8), (10). | 



5 Conclusion and Acknowledgements 

Univariate polynomials with nonnegative real roots appear quite often in modern combinatorics 
, especially in the context of integer polytopes . The closest to our approach is the class of 
univariate rook polynomials [2HI • We discovered in this paper rather unexpected and very 
likely far-reaching connections between hyperbolic multivariate polynomials and many classical 
combinatorial and algorithmic problems . 

The main "spring" of our approach is that the class of PCS'- hyperbolic polynomials is large 
enough to allow the easy induction . The reader might be surprised by the absence of Alexandrov- 
Fenchel inequalities and other ingredients of proofs in , j i31j . In fact , the clearest (in 
our opinion) proof of the Alexandrov-Fenchel inequalities for mixed discriminants is in A.G. 
Khovanskii' 1984 paper • The Khovanskii' proof is based on the similar induction (via 
partial differentiations) to the one used in this paper . In a way , the Alexandrov-Fenchel 
inequalities are "hidden" in our proof . 
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Let us summarize the main ideas of our approach : 

IDEA 1 . To facilitate the induction we deal not with doubly stochastic matrices/tuples/polynomials 
but rather with the CAPACITY of homogeneous polynomials with nonnegative coefficients . 
IDEA 2 . The notion of the VDW-FAMILY allowed to reformulate Van der Waerden / 
Schrijver-Valiant/Bapat conjectures in terms of homogeneous polynomials with nonnegative 

coefficients . 

IDEA 3 . The notion from the theory of linear PDE , POS'-hyperbolic polynomials , hap- 
pened to give the needed , i.e. containing multilinear and determinantal polynomials , VDW- 
FAMILY . Corollary 2.8 , a particularly easy case of the Van der Waerden Conjecture , was 
the final strike . 

I would like to thank Sergey Fomin , Shmuel Friedland , Mihai Putinar, Alex Samorodnitsky 
, Lex Schrijver , Warren Smith , George Soules , Ian Wanless , Hugo Woerdeman for the interest 
to this paper . 

This paper was , to a great degree , Internet assisted (sometimes in mysterious ways). Thanks 
to Google ! 
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1. Doing simple "algebra" we get that 



Sn-l-nSn= Ci( ^ —). 

l<i<n l<i<n 

Notice that < 1 — q < 1 and J2i<i<ni^ ~ ^i) — 1- Using the concavity of the logarithm 
we get that 

logiSn-i -nSn)> E log(ci)+ E (l-Q)log(-)= ^ (^Md). (20) 

l<i<n l<i<n ^ l<i<n 

2. 

l<j<n j^i 

Define Cj = 1 — Oj . Then < 1 — Cj < 1 , Z^i<j<„ Cj = ra — 1 and the permanent 

= (i - l)n-l (^''-l - 

It is easy to prove and well known that 

min Cilog{ci)= E -^og{- i-) = log((^^^^ 

Using the entropic inequality (20) from the first part we get the following equality 

min Sn-i-nSn = {'^r-'. 

0<1-Ci<l;^j^.^„c,=n-1 n 

Which gives the needed inequality 

PeriA) > j^^^^i^r-' = -. (21) 

It is easy to see (strict concavity of Z]i<i<„ Q log(ci)) that the last inequality is strict 
unless A{i,j) = ^; 1 < i, j < n. 



B Proof of Corollary 2.8 

Proof: We can assume WLOG that Qi > 0,1 < i < n.If R{0) = then clearly di = R'{0) > 

C > C((^^)"-i) . Let R{0) > , i.e. bi> 0,1 <i <n. 

Associate with polynomial R{t) = ni<j<n(^i* + ^i) the following matrix with positive entries 
A = [a\c\...\c] , where a = (ai, a^)^, c = ;^(6i, bn)'^- The condition R(t) > Ct,\/t > is 



17 



equivalent to the inequality Cap{MulA) > C . And per{A) = j^z^^^di. Since A has positive 
entries hence there exist two positive diagonal matrices Di,D2 such that A = D1BD2 , where 
the matrix B is doubly stochastic and B = [f\d\...\d] (Sinkhorn Diagonal Scaling , see Lemma 3.6 
in ^] in a more general setting). Since the matrix B is doubly stochastic hence Cap{MulB) = 1 
. Thus Cap{MulA) = det{DiD2)Cap{MulB) = det(L>iL>2) and per{A) = det{DiD2)per{B). 
Therefore , we get from inequality (21) that per{A) > ^Cap{MulA) > ^C. 

Finally ,it follows that 

[n — Ij" ^ (n — Ij" ^ n" n 



C Proof of Theorem 2.5 

We need the following simple result . 

Proposition C.l: Let p{X) be e-hyperbolic (homogeneous) polynomial of degree n , p{e) > 

. Consider two e-nonnegative vectors Z,Y N(,{p) such that Z + Y ^ Ce{p) , i.e. Z + Y is 
e-positive . Then 

p{tZ + Y)= n M + > 0,ai + 6i > 0,1 < i < n. (22) 

l<i<n 

Proof: As the vector Z + Y = D \s e-positive hence p(Z-|-y) > (FACT 1) ,the polynomial p 
is Z + y- hyperbolic (FACT 3) and any e-positive (e-nonnegative) is also Z-|-y-positive(Z-|-y- 
nonnegative) (FACT 3) . Doing simple algebra , we get that p{tZ + Y) = p{{t — 1)Z + D). 

Let < Ai < A2 < ... < An be nonnegative roots of the equation p[Z — xD) = 0. Since 
D - X = Y ^ Ne{p) = Nd{p) hence A„ < 1. Therefore 

p{tZ + Y)=p{{t-l)Z + D)=p{D) II (tA, + (l-A,)) 

l<i<n 

We can put = (p(Z + Y))\i >0,bi = {p{Z + Y)){1 - Ai) > 0. I 

C.l Proof of the first part of Theorem 2.5 - inequahty (16) 

Proof: Let q € Homj^{n, n) be P05-hyperbolic polynomial and 1 < Rankq{ei) = Sq{{l}) = k. 
Fix positive real numbers {x2, x^) such that Yl2<i<n^i ~ 1- Define the following two real 
n-dimensional vectors with nonnegative coordinates : Z = (1, 0, 0, 0), F = (0, X2, x„) . 
The vector Z -|- y is e-positive . Consider the next univariate polynomial R{t) = q{tZ + Y). It 
follows from Proposition C.l that 

m= n M+&i)= E 

l<i<n 0<i<n 
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where ai,bi > and qxi{x2, = di and the cardinahty \{i : ai > 0}\ = k (see also equality 

(2)). In other words the degree deg{R) = k 

We get straight from the definition of Cap{q) that 

Rit) = Yi {o-it + bi) =p{t,X2,...,Xn) >Cap{q)t J| Xi = tCap{q). 

l<i<n 2<i<n 

Using Corollary 2.8 , we get that 

A" — 1 J. 

qxi{x2,--,Xn) =di> ((— ^) ~^)Cap{q). 
In other words , that Cap{r) > {^)'^~^)Cap{q). I 

C.2 Proof of the second part 

The second part of Theorem 2.5 is an easy modification of FACT 4. We need only to consider 
the 7^ 0. We need the following well known fact . 

Fact C.2: Consider a sequence of univariate polynomials of the same degree n : Pkit) = 
Eo<i<n (^iM'' ■ suppose that limfc^oo ai,fc = ai,0 <i <n and 7^ . 

Define P{t) = J2o<i<n ^i^^ ■ Then roots of Pfc converge to roots of P . In particular if roots of 
all polynomials are real then also roots of P are real ; if roots of all polynomials P^ are real 
nonnegative numbers then also roots of P are real nonnegative numbers . 

It follows from Definition 2.2 and the Taylor's formula (14) that 

g<i(X) = ^g(X + td)|(t=o) = ((^-l)!)-'M,(d,X,X,...,X) :d,XGP" (23) 

Notice that qx^{x2, ■■■■,Xn) = ^ei (0, 2^2; a^n) ■ Consider the following perturbed univariate 
polynomials P^it) = Eo<j<n-i 

Pe{t) = ((n - ly.y^Mgiei + ee,Y- t{{e - d) + ee), Y - t{{e - ei) + ee)) : 
e= J2 ei,y = (X2,...,x„) G P"~\e > 0. 

l<i<n 

We get by a direct inspection that 

^lim Pe(t) = ((n - iy.)-^Mqiei,Y - t{e - d), Y - t{e - ei)) = q^^ (xa - t, x„ - t); 

lim ae,n-i = (-l)"~^fe(l, 1, !)• 

As qxi G M_|_(n — 1, n — 1) and qx^ 7^ hence qxi{^, 1, 1) > . Therefore lim£_>oo ae,n-i 7^ . 
Since the polynomial q G Hom+{n,n) is POS"- hyperbolic hence it follows from FACT 3 that 
all the roots of the equation Pe(t) = are real ; if > 0, 2 < i < n then all the roots are 
nonnegative . We conclude, using Fact C.2 , that 
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1. If y = (x2, ■■■,Xn) G ^ then all the roots of the equation QxiiY — t{e — ei)) = are 
real and qxi{e — ei) > 0. 

2. If Xj > 0, 2 < i < n then all the roots are nonnegative . 

3. Since € M+(n — l,n — 1) and ^ hence qxi{x2, Xn) > if > 0, 2 < z < n . 
It follows from the equality (15) that all the roots of the equation QxiiX — t{e — ei)) = 
are positive if Xj > 0, 2 < i < n . 

4. The polynomial qx^ G M+(n — l,n — 1) is P05-hyperbolic . 

D Proof of Theorem 4.6 

Proof: Let q € Hom^{n,n) be P05-hyperbolic polynomial and a pair of indeces (i ^ j) C 
{1, 2, n} . Define the following integer vectors A^'^^ = e + — ej. Condition 1 states that 
j-i^J) £ supp{q) for all such pairs . The equivalence of Condition 1 and Condition 2 follows 
from the second part of Theorem 4.1. 

The fact that Condition 1 implies indecomposability is valid for all polynomials in Homj^[n, n) 
and is proved in pLlj] . 

Suppose that there exists a positive vector /? = (3n),I\i<i<n Pi — ^ such that 

q{(]i,...,(3n) = inf q{xi, Xn) = Cap{q). 

Then the polynomial Q{xi, ...,Xn)) = ^ /; ) Q{J2i<i<n Pi^i) is doubly-stochastic . Notice 
that the inf^ . q tt qixi, ■■■,Xn) is attained and unique if and only if the inf^.. q tt ^^^^ 

is attained and unique ; it follows from the Euler's identity that inf^.^Q tt ^.^-^ Q{xi^ ^n) = 
1 = Q{e) . _ 

Assume that Condition 2 does not hold : there exists a subset A C {1, 2, n}, 1 < = m < 
n such that Rankg(J2ieA (^i) = R'in'kQ(J2ieA(^i) = ^- Define ca = T^i^^A^^ii^^A' = Y^jeA'^^j = 
e — SA- Let Xn > ■■■ > A„_m+i > > ... > be the ordered roots of the equation Q{eA — te) = 
. Since ba = e — ca' and the vector ca' is e-nonnegative hence 0<Ai<l,l<i<n. As the 
polynomial Q G Homj^{n,n) is doubly-stochastic hence (see ^H] i E3) 

^ A, = ^ X-=m= \A\. 

l<.i<n n—m+l<i<n 

Therefore Xi = l,n — ni + l<i<n and Xj = 0,1 < j < n — m . It follows from the identity 
(15) that Q{aeA + beA') = 1 = Q{e) iff a™-5"-™- = l ^ which proves the non-uniqueness of 

The following result is proved very similarly to the previous proof , the only new ingredient 
is the subadditivity of Rankq{X), X E R". It "justifies" the notion of " indecomposibility of 
POS'-hyperbolic polynomials" . 
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Theorem D.l: Let q G Hom+{n,n) be POS -hyperbolic polynomial . Supposed that 

inf^^^Q TT ^.^^ q{xi, Xn) is attained and Rankg{J2i<i<m ^i) = m,m < n. 
Then the polynomial q can be decomposed in the following way : 

= qiixi, ..,Xm)q2{xm+i, -.Xn) • qi e Hom+{m,m),q2 G Hom+{n-m,n-m) 

(24) 
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